Doping dependence of the Fermi surface in (Bi,Pb)2Sr2CaCu208+5 
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A detailed and systematic ARPES investigation of the doping-dependence of tiie normal state 
Fermi surface (FS) of modulation-free (Pb,Bi)-2212 is presented. The FS does not change in topology 
away from hole-like at any stage. The FS area does not follow the usual curve describing Tc vs x 
for the hole doped cuprates, but is down-shifted in doping by ca. 0.05 holes per Cu site, indicating 
the consequences of a significant bilayer splitting of the FS across the whole doping range. The 
strong fc-dependence of the FS width is shown to be doping independent. The relative strength of 
the shadow FS has a doping dependence mirroring that of Tc. 

PACS numbers: 74.25.Jb, 74.72. Hs, 79.60.-i, 71.18.+y 



The shape and topology of the Fermi surface (FS) 
of the high temperature superconductors (HTSC), and 
in particular of the Bi2Sr2CaCu208+5 (Bi-2212)-based 
systems, has been a hat-, topic from the very begin- 



ning of the HTSC era 
lively discussion today.[ 



and is still the subject of 
In the past, the existence 
of a large, hole-like FS in angle-resolved photoemission 
spectroscopy (ARPES) was taken as support for the va- 
lidity of jJjittinger's theorem for the superconducting 
cupratesJa-El While some ARPES studies of Bi-2212 con- 
clude that a Urge, hole-like FS persists even to very lo- 
doping levels,El other data imply a change in FS topolog; 
or the presence of hole-pockets at underdoping.E2l Recent 
data from La2-2;Sra;Cu04 (LSCO) have been interpreted 
in terms of a change of FS topology from hole-like for 
X < 0.2 to electron-like for higher doping levels.tll 

The recent improvement in the performance of pho- 
toemission instrumentation (in particular in the angular 
resolution) has led to a renaissance in the direct deter- 
mination of the basal plane projection of the FS using 
ARPES. Considering the fundamental importance of the 
FS topology and shape in deciding the physical proper- 
ties of a solid, it is natural to want to study its doping 
dependence in the Bi-based HTSC directly and with high 
precision using high resolution FS mapping. 

The ARPES experiments reported here were per- 
formed either using monochromated He I radiation and 
an SES200 electron analyzer or using hv = 25 ey-radi- 
ation from the U125/1-PGM beamhne at BESSYB and 
an SESlOO electron analyser. The samples were mounted 
on a triple-axis goniometer, enabling computer controlled 
angular scanning with a precision exceeding 0.1° for all 
axes, resulting in a dense sampling of a large portion 
of kcj-space for each single crystal studied. The overall 
resolution was set to 0.014 A~^x 0.035 A'^x 19 meV 
which are the FWHM momentum (parallel and perpen- 
dicular to the analjper entrance slit) and energy reso- 
lutions, respectively.li3 The samples were cleaved in-situ 



to give mirror-like surfaces and all data were measured 
above the pseudogapped regime at 300K within 3-4 hours 
of cleavage. The synchrotron based data were collected at 
30 K with 0.014 A'^x 0.014 A'^x 17 meV resolution. 
We investigated a set of high quality single crystals of 
Pb-doped Bi-2212 (PbrBi ratio = 0.4:1.6) which had un- 
dergone different oxygen loading procedures. As we have 
pointed out earher (see, e.g., Rcfs. p| Jl4|p^ ), it is wise to 
use the Pb-substituted variants for such experiments as 
these systems do not possess the incommensurate modu- 
lation of the BiO layers which in pristine Bi-2212 leads to 
the appearance of strong diffraction replicas of the main 
and shadow FS features in the maps, thus disqualifying 
a detailed discussion of the FS topology, shape and area 
as a function of doping. In the following, we label the 
samples, which span a range of 35 K around optimal 
doping, according to their T^. UD 76K, UD 85K, UD 
89K, CD 81K, OD 72K and CD 69K (UD and OD stand 
for undcrdoped and overdoped). 

Fig. |l| shows the Fermi surface maps for all six doping 
levels. Each dataset contains ca. 5000 ARPES spectra. 
We collect data from a significantly larger region of k- 
space than the irreducible octant, which brings the ad- 
vantage of enabling a quantitative correction of angular 
misalignments of the crystal to a precision of 0.1°. 

To minimize the effects of the factors separating the 
ARPES intensity distribution from the spectral function 
(see Ref. |l3|), the data were 'self- normalized' by divid- 
ing the signal from the Fermi level, /(k,Lj = 0), by the 
signal at highest binding energy, /(k, w/ifce) (here uJube 
is 300 meV). The FS topology and shape derived from 
these data do not depend sensitively upon the use of any 
reasonable self-normalization denominators, although we 
wish to stress here that the self-normalization procedure 
itself is indispensable for the precise determination of the 
kf-vectors (see Appendix). 

Before going on to discuss the data in a more quantita- 
tive manner, we first cover what can be learned directly 
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FIG. 1: Basal plane projection of the normal state (300K) Fermi surface of Bi(Pb)-2212 from high resolution ARPES. The Ef 
intensity (normalized to the signal at a; = 0.3 eV) is shown in color. The Tc of each sample is indicated. The raw data cover 
half of the coloured area of each map, and have been rotated by 180° around the F point to give a better k-space overview. 
The line dividing raw and rotated data runs almost vertically for the UD76K map, and from top left to bottom right in all 
other maps. The sketch shows the FS for the OD 69K dataset as yellow barrel-like shapes defined by joining the maxima of fits 
to the normalized Ef MDC's. 



from a simple visual inspection of Fig. |lj. (i) There is 
no topological change of the main FS within the dop- 
ing range studied - it remains hole-like (centered at the 
X, Y (±7r,±7r) pOjiiLts), in contrast to recent data from 
the LSCO systemJld (ii) As hole doping is increased, the 
main FS 'barrels' increase in size (as can easily be seen 
in the decrease of the inter-barrel separation around the 
M (tTjO) point), accompanied by an increase in the size of 
the lenses formed by main FS and shadow FS (SFS). (iii) 
The shape of the FS barrels changes from being quite 
rounded at low doping to take on the form of a square 
with well-rounded corners at higher doping, (iv) The SFS 
exists at all doping levels. 

We stress that these statements describe experimental 
observations and are independent of any particular data 
analysis or physical interpretation. 

One of the fundamental questions in the physics of 
2D strongly correlated electron systems is to what ex- 
tent the interacting electron system can be described by 
models derived perturbatively from the non-interacting 
case. One way to test this is to consider the validity 
or otherwise of Luttinger's theorem, which can be para- 
phrased by stating that the volume (area in 2D) of the FS 
should be conserved upon switching on the interactions. 
Thus if we are able to pin down the doping dependence 
of the exact path in k-space which represents the Fermi 
surface in, for example, the (Pb,Bi)-2212 HTSC without 
knowing, a priori, its shape, we would be able to evalu- 



ate the doping dependence of the FS area and thus test 
Luttinger's theorem. The best approach here is to locate 
the maxima in the Ep momentum distribution curveslij 
(MDC's) describing tracks crossing the FS (preferably 
at right angles, see Appendices here and in Ref. |l^ for 
details). 

Such a fitting procedure was carried out for the OD 
69K sample. The detailed result is well described by a FS 
having the form of a square with rounded corners, which 
confirms the visual impression from the intensity map 
for this sample. Such a form gives a simple analytical 
approximation for the FS shap e als o predicted in tight- 
binding and LDA calculations .EZlta A sketch of the fit 
result is shown as the yellow line on the right hand side 
of Fig. |. The FS maps from the other samples were 
then fitted, whereby the extent of the straight sections, 
as well as the size of the barrel as a whole were varied to 
optimize the fit to the data. We can then derive the hole 
concentration x from the simple relation x+1 = 2Sb/ Sbz, 
where Sb is the area of main FS barrel and Sbz is the area 
of the Brillouin zone. 

The results obtained from the analysis of the FS area 
are shown in Fig. |^ in the form of a Tc vs a; plot. The solid 
line shows the commonly employed empirical relation be- 
tween Tc and xS3 For the six samples spanning a total 
of 35K in Tc, the co-ordinate pairs matching the T^s to 
the doping level taken directly from the experimentally 
determined Fermi surface area also give a parabolic curve 
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FIG. 2: Symbols: critical temperatures vs the hole con- 
centration, xfs, the latter beinecalculated directly from the 
area of the FS's shown in Fig. |l]. The solid line shows the 
commonly-used empirical relation for Tc vs x (Ref. |l|). 

(shown as a dotted line), but this curve is down-shifted 
in doping by ca. 0.05 towards the underdoped side of the 
phase diagram. This result, being quite surprising not 
long ago, can be well understoodjpw in terms of the 
bilayer splitting of the CuO bandES-El 

Before going further, we note that such a shift is hard 
to explain by the assumption that the doping level at 
the surface is lower than in the bulk. If this were the 
case (for example by loss of oxygen at the surface), such 
a deviation should be strongly dependent on the oxy- 
gen loading procedure, affecting the OD samples more 
strongly than the UD, which is clearly not the case. Fur- 
thermore, the fact that the superconducting gap seen in 
ARPES data from the same samples (not shown), closes 
unambiguously at the bulk Tc in the overdoped systems, 
is incompatible with a lower doping level at the surface. 

The acceptance of an existence of the c-axis bilayer 
splitting in Bi-2212 marks a watershed in the interpreta- 
tion of ARPES data from the multilayer HTSC. This 
splitting has beenj-dkectly resolved recently in highly 
overdoped Bi-2212BEl and (Bi,Pb)-2212El and shown 
(but notxesolved) to be roughly the same for underdoped 
Bi-2212B. In Fig. H, we show azimuthal energy distribu- 
tion maps (EDM's: I{d,uj), where 6 is the azimuth angle, 
see Fig. ||) at T = 30 K for three different |k| = 1.084, 
1.088 and 1.092 A^^ (from the left to the right corre- 
spondingly) which demonstrate a well resolved bilayer 
splitting in underdoped (Bi,Pb)-2212 {Tc = 77 K). 

Given the presence of the bilayer splitting (which we 
include here in the notion 'complex structure'), the blur- 
ring of the FS (see Fig. |l|) on going from the nodal to 
the antinodal point for all doping levels, which is often 
attributed to the complex physics of antinodal electrons 
(e.g. an absence of well-defined quasiparticles) , could at 
least partially be due to such a complex structure of the 
FS itself. In order to examine this possibility, in the fol- 
lowing we analyze the FS width in more detail. 
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FIG. 3: The azimuthal energy distribution maps (EDM's) at 
T = 30 K for three different |k| = 1.084, 1.088 and 1.092 A"^ 
(from the left to the right correspondingly) which demonstrate 
a well resolved bilayer splitting in underdoped (Bi,Pb)-2212 
in = 77 K). 

In Fig. m we show the width of the FS, Ak vs 0, the 
latter being the angle away from the nodal line, as indi- 
cated in Fig. |. The Afc values were derived from fitting 
Ep MDC's using a Lorentzian profile with Ak FWHM. 
For all doping levels investigated the room temperature 
FS width is strongly fc-dependent, being maximal near 
the antinode and minimal at the node. The dotted line 
in Fig. 1^ shows that the data can be well described by 
the function 

A/c((?!)) = A/co + Afcisin2(2(/)), (1) 

where Afco = 0.054 k'^ and Afci = 0.136 A^^. 

Remarkably, the observed fc-dependence of the FS 
width is essentially independent of the doping level. This 
is difficult to reconcile with a FS width determined solely 
by the complex physics of the FS electrons, as within 
such a picture the difference in the coupling to inter- 
actions between the nodal and antinodal regions should 
decrease continually as the doping increases. Equally, we 
can rule out effects resulting from differing group veloci- 
ties around the FS contour, as these have been shown to 
be essentially constant .E3 

On the other hand, exploiting the 'complex FS struc- 
ture' scenario, we can associate a splitting in momentum, 
dk{(j)), with the c-axis bilayer splitting. For the case in 
which the maxima of the MDCs (i.e. the intensity in a 
self-normalized FS map such as those of Fig. |l|) corre- 
spond to the inner bilayer split FS barrel (namely the 
bonding CuO-bilayer band), this would result in a shift 
of the observed doping level of 

r 5km<l^, (2) 

OBZ JBZ Jo 

where 5S is the difference in area between the split bar- 
rels, kb is the radius of main FS barrel with respect to 
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FIG. 4: The width of the room temperature main FS, AA: vs 
the FS angle (j> defined with respect to the nodal line. The Tc's 
are indicated and the solid grey line represents the relation 
(0), for details see text. 



the X-point ((fcb) « 0.6|rX|) and |rX| = 1.161 A^^. This 
effect is illustrated schematically in the cartoon shown in 
Fig. where the yellow (red) barrels represent the smaller 
(larger) FS's resulting from the bilayer splitting. Taking 
a Lorentzian form for the Ep MDC which cuts the FS, 
we expect a bilayer splitting induced FS width given by 



Ak^W 



3{Skf 
2W 



(3) 



where W is the FWHM of the FS without splitting and 
5k < W/VS is assumed to hold). In such a manner we 
can estimate an upper limit for Sx — 0.07, which is illus- 
trated in Fig. H by the broad grey arrow. This demon- 
strates that the effect of the bilayer splitting is enough 
to explain the downshift of the vs doping parabola. 

The reason why the photocurrent intensity from the 
antibonding band at Ep is less than from the bonding 
one and, consequently why the maxima of the MDCs 
correspond to the inner bilayer split FS barrel is the dif- 
ference between matrix elements for photoemission from 
these two bands. In fact, the ratio between the effec- 
tive matrix elements for emission from the bonding and 
antibonding Cu-0 bands, Ah/Ma « 2 for hv=21.2 eV 
(see Fig. 3 in Ref. |l^). In contrast, for 25 eV excitation 
energy, Mb/ Ma ~ 1 and neither the emission from the 
bonding nor the antibonding band dominates resulting 
in a clear splitting as can be seen in the EDM's shown in 
Fig.|. 

Finally, we note in this context that the upper limits of 
5k (and consequently 5x) obtained above correspond to 
the limit at which two Lorentzian features are resolvable 
from one another: 5k = W/^/i. This same limit also de- 
fines the lower bound for the 0-dependence of the Fermi 
surface width which arises from sources other than the 
bilayer splitting: 



Afc(0) = Afco(l + 1.3sin2(20)). 



(4) 



In other words, this means that the detected anisotropy 
described by Eq. ([|) cannot be explained by the bi- 
layer splitting alone. In considering either the 'complex 
physics' or 'complex FS structure' scenarios we discuss 
two extremes, whereas the real situation may well in- 
clude contributions from both. For example, at high hole 
doping, the ^-dependence of Afc from 'complex physics' 
should flatten out, which would be counteracted by the 
increasing bilayer splitting for this doping regime (in 
which the flat bands approach closer to Ep)- Conversely, 
at low hole doping, the 0-dependence of the coupling to 
interactions is strong, wheiieas the bilayer splitting would 
be expected to be weaker .tZI In this way we end up with 
the observed overall doping independence of Ak{(j)). 

As mentioned above, it is possible to compensate for 
the downshift of the Tc vs. x parabola in Fig. 2 by tak- 
ing the bilayer splitting into account. It would then fol- 
low that the area of the main ARPES FS scales with 
(1-l-a:) in holes across the complete doping range studied. 
This behaviour is in contrast to what is seen in transport 
measurements. Resistivity and Hall effect datsuipdicate 
that the transport characteristics scale with x,E3'E3 even 
into the overdoped regimeEj Although it is conceivable 
that only those mobile electrons which have relatively 
low coupling to other degrees of freedom contribute to 
the transport, it is surely more than coincidental that 
this proportion should be exactly x/{l+x). This fun- 
damental difference between the transport data and the 
ARPES FS is a key question which deserves detailed the- 
oretical attention. 

A final surprise that the FS has in store for us is shown 
in Fig. 1^ (lower panel), in which the doping dependence 
of the intensity ratio of the SFS to that of the main FS 
is plotted. The intensities Im were taken in each case 
from the same azimuthal MDC scan: i.e. with the same 
|k| value, some 0.13 from the point at which the 
SFS and main FS 'cross'. For this scan the given in- 
tensity ratio reaches a local maximum as a function of 
|k| which is a consequence of different dependencies of 
the photocurrent from the main and shadow FS's (at 
the same |k|) on matrix elements. The upper panel of 
Fig. 1^ shows an example (for the OD 69K) of such az- 
imuth MDCs from which these intensity ratios have been 
determined: peaks 1 and 4 correspond to the SFS, peaks 
2 and 3 correspond to the main FS, then the intensity 
ratio SFS/FS = h)/{l2+ h)- As Fig. | shows, this 
ratio decreases not only on going from optimal to over- 
doping, but also on going towards the underdoped side 
of the phase diagram (the rate of change is, in fact, even 
faster on the UD side). This is in contrast to predictions 
based on an antiferromagnetic origin of the SFS. El The 
fact that the relative strength of the SFS tracks the dop- 
ing dependence of Tc means that, regardless of whether 
the SFS has structural or other origins, this phenomenon 
is important and could be related to high Tc supercon- 
ductivity itself. The behaviour seen here, taken together 
with the very strong similarities with SFS data from pris- 
tine Bi2212 (which has important structural differences 
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combination of the effects of the bilayer spHtting (domi- 
nating at higher doping) and the complex physics of the 
FS electrons (dominating at lower doping). Finally, the 
shadow FS is clearly visible for all doping levels, and has 
maximal intensity at optimal doping, raising the ques- 
tion of a possible link between the origins of the shadow 
FS and superconductivity. 
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O. K. Andersen for stimulating discussions. 
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FIG. 5: SFS to main FS intensity ratios vs T^"'' - Tc (lower 
panel), the dashed straight lines are guides to the eye; and 
an example of an azhimuth MDC from which these intensity 
ratios have been determined (upper panel): peaks 1 and 4 
correspond to the SFS, peaks 2 and 3 correspond to the main 
FS, then the intensity ratio SFS/FS = [h + h)/{h + h). 



to Pb-doped Bi2212) means that further work is needed, 
both on the experimental but also on the theoretical side, 
before the question of the origin and consequences of the 
shadow Fermi surface can be considered as being solved. 

In conclusion, we have presented a detailed and sys- 
tematic ARPES investigation of the doping dependence 
of the normal state (room temperature) FS of the Bi- 
2212 family of HTSC materials. The data clearly show 
no change in the FS topology away from hole-like at 
any stage (from UD 76K to OD 69K). An analysis of 
the main FS area gives a parabolic Tc vs xps relation, 
shifted to lower x by some 0.05 compared to the 'univer- 
sal' relationjlj which can be accounted for by the pres- 
ence of two (unresolved) FS's near (7r,0) due to a bilayer 
splitting with a maximum value ca. 0.05 A~^, which stays 
roughly constant across the whole doping range. Further- 
more, the FS width is shown to be strongly dependent on 
k, but for each particular kp point it is essentially inde- 
pendent of the doping level. This can be understood as a 



APPENDIX A 

A meaningful estimation of the doping level from a FS 
map as described above requires the very precise determi- 
nation of the FS vectors. Here we explain why the self- 
normalization procedure has been chosen for this pur- 
pose. 

In Ref. we already discussed the applicability of dif- 
ferent methods of fc^ determination to ARPES data from 
Bi-cuprates and demonstrated-that the most accurate is 
the 'maximum MDC method. Ej In Ref. |l^ this was illus- 
trated for the case of the nodal direction. At lower energy 
resolution (or for the case in which the MDC peaks are 
broader), however, the deviation of the experimentally 
determined kp from the true value, Afci?, could be consid- 
erable, and even comparable with the A/c^^'s from other 
methods such as 'gradient n(fc)' (see Ref. ^3|). It turns 
out that the above mentioned shift (Afc^?) is nearly com- 
pletely compensated for by the self-normalization proce- 
dure and therefore the deviation of the visible FS traces 
on the self-normalized intensity maps (like those shown 
in Fig. 1^) from the real FS is negligible. This is demon- 
strated below. 

Fig. H represents the results of Ep-MDC simulations in 
the nodal (F— (tt, tt) crossing, upper panel) and antinodal 
(((tt, 0) — (7r,7r) crossing, lower panel) points for typical 
low-energy dispersion relations — vp{kp — k) with 
vp = 2 eVA. For this simulation we use a simple form 
for the spectral function (with the momentum resolution 
included) 



A{k,U!, Rk) cx 



v/S"(c.,T)2 



(w-efc)2 + E"(cj,T)2 + 
where the imaginary part of the self-energy 
S"(tj,T) 



Vic 



{pry 



(Al) 



(A2) 



has been taken with a = 1 and (3 = 2 (binding energy w 
and temperature T = 300K are i n l en ergy units), giving 
a resonable fit to the experimentOEJ. The real part of 
the self-energy is included in e^. The photocurrent is 
calculated as 



I{k,Lo) cc [A{k,u,Rk)f{L^)](g)R^, 



(A3) 
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where /(w) is the Fermi function. The momentum res- 
olutions are taken to be 0.014 A^^ and 0.035 A^^ for 
the rX and MX crossings respectively and energy res- 
olution of 19 meV for both. For the MX crossing, to 
check an extreme, the 'splitting value' Ae — 80 meV is 
added to the R^^ function as a FWHM. The dashed curves 
represent the raw (non-normalized) MDCs: /(fc,0). The 
solid curves represent MDCs obtained by self-normalizing 
every EDC to the highest binding energy {ujhbe in this 
case is 0.3 eV): I{k,0)/I{k,LUiihi.)- For comparison, the 
dn{k)/dk dependencies are shown as dotted lines where 

I{k,uj)dLu. (A4) 

Fig. ^ illustrates that whereas in the case of the nodal 
region both raw and self-normalized MDCs are only 
slightly shifted from the real kp {Akp — —0.001 and 
0.003 A~^ respectively), in the antinodal region the shift 
of the raw MDC is rather large (—0.015 A~^) whereas 
the peak of the self-normalized curve practically coincides 
with the true kp (i.e. Akp = -0.001 A^^). This demon- 



strates the power of the self-normalization procedure: its 
application to the intensity maps not^only reduces the 
influence of the matrix element effectala but also restores 
the true location of the FS vectors, thus making it the 
correct choice in the study of the FS topology, shape and 
area. 
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FIG. 6: Results of Ef-MDC simulations for nodal (F - (tt, tt) 
crossing, upper panel) and antinodal ((tt, 0), (tt, tt) crossing, 
lower panel) ARPES data for a typical low-energy dispersion 
(for details see text). The dashed curves are the raw (non- 
normalized) MDCs, the solid curves are the MDCs after self- 
normalization to the highest binding energy, and the dotted 
curves represent the results of the dn{k)/dk method. 
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